In vitro analysis of Porcine Endogenous Retroviruses in different pig cell types.

The shortage of organs for human transplantation is a topic of extreme interest, and xenotransplantation with porcine organs has been recognized as a promising solution. However, the potential spillover linked to infectious agents present in pigs remains a concern. Among these, Pig Endogenous Retroviruses (PERVs), whose proviral DNAs are integrated in the genome of all pig breeds, represent an extremely important biological risk. This study aims to evaluate PERVs distribution in several swine cell lines and samples of domestic and feral pigs. Moreover, the capacity of PERVs to infect human and non-human primate cells and to integrate in the cellular genome was tested by Real-Time PCR and by Reverse Transcriptase assay. Results indicated a widespread diffusion of PERVs both in cell lines and samples analysed: the viral genome was found in all the established cell lines, in 40% of the primary cell lines and in 60% of the tissue samples tested. The assays indicated that the virus can be transmitted from porcine to human cells: in the specific case, infected NSK and NPTr cells allow passage to human 293 and MRC-5 cells with active production of the virus demonstrable via PCR and RT assay. In light of these aspects and also the lack of studies on PERVs, it appears clear that there are still many questions to be clarified, also by means of future studies, before xenotransplantation can be considered microbiologically safe.

A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record.

BACKGROUND: Genomic regions that remain poorly understood, often referred to as the dark genome, contain a variety of functionally relevant and biologically informative features. These include endogenous viral elements (EVEs)-virus-derived sequences that can dramatically impact host biology and serve as a virus fossil record. In this study, we introduce a database-integrated genome screening (DIGS) approach to investigate the dark genome in silico, focusing on EVEs found within vertebrate genomes. RESULTS: Using DIGS on 874 vertebrate genomes, we uncover approximately 1.1 million EVE sequences, with over 99% originating from endogenous retroviruses or transposable elements that contain EVE DNA. We show that the remaining 6038 sequences represent over a thousand distinct horizontal gene transfer events across 10 virus families, including some that have not previously been reported as EVEs. We explore the genomic and phylogenetic characteristics of non-retroviral EVEs and determine their rates of acquisition during vertebrate evolution. Our study uncovers novel virus diversity, broadens knowledge of virus distribution among vertebrate hosts, and provides new insights into the ecology and evolution of vertebrate viruses. CONCLUSIONS: We comprehensively catalog and analyze EVEs within 874 vertebrate genomes, shedding light on the distribution, diversity, and long-term evolution of viruses and reveal their extensive impact on vertebrate genome evolution. Our results demonstrate the power of linking a relational database management system to a similarity search-based screening pipeline for in silico exploration of the dark genome.

HERV-W upregulation expression in bipolar disorder and schizophrenia: unraveling potential links to systemic immune/inflammation status.

BACKGROUND: Bipolar disorder (BD) and schizophrenia (SZ) are the two main mental disorders with unknown etiology that significantly impact individuals' quality of life. The potential pro-inflammatory role in their pathogenesis is postulated and Human Endogenous Retrovirus W (HERV-W) is an emerging candidate to modulate this pathogenic finding. HERVs, ancient retroviruses in the human genome, may play roles in inflammation and disease pathogenesis. Despite HERVs' involvement in autoimmune diseases, their influence on mental disorders remains underexplored. Therefore, the aim of this study was to assess the level of HERV-W-env expression and the systemic inflammatory profile through the concentration of IL-2, IL-4, IL-6, IL-10, TNF-α and INF-γ cytokines in BD and SZ patients. RESULTS: All participants showed HERV-W-env expression, but its expression was higher in mental disorder patients (p < 0.01) than in control. When separated, SZ individuals exhibited higher HERV-W expression than the control group (p < 0.01). Higher serum levels of TNF-α and IL-10 were found in BD (p = 0.0001 and p = 0.001, respectively) and SZ (p = 0.01) and p = 0.01, respectively) than in the control group, while SZ showed decreased levels IFN-γ and IL-2 as compared to controls (p = 0.05) and BD patients (p = 0.05), respectively. Higher TNF-α/IL-4 and TNF-α/IL-10 ratios, and lower IFN-γ/IL-10 were observed in BD and SZ patients than controls. Significant negative correlation between HERV-W-env expression and IL-10 (r=-0.47 p < 0.05), as well as positive correlations between HERV-W-env expression and TNF-α/IL-10 or IFN-γ/IL-10 ratios (r = 0.48 p < 0.05 and r = 0.46 p < 0.05, respectively) were found in BD patients. CONCLUSION: These findings suggest not only a potential link between HERV-W-env expression both in BD and SZ, but also a possible involvement of systemic inflammatory status in BD patients.

Adaptive expansion of ERVK solo-LTRs is associated with Passeriformes speciation events.

Endogenous retroviruses (ERVs) are ancient retroviral remnants integrated in host genomes, and commonly deleted through unequal homologous recombination, leaving solitary long terminal repeats (solo-LTRs). This study, analysing the genomes of 362 bird species and their reptilian and mammalian outgroups, reveals an unusually higher level of solo-LTRs formation in birds, indicating evolutionary forces might have purged ERVs during evolution. Strikingly in the order Passeriformes, and especially the parvorder Passerida, endogenous retrovirus K (ERVK) solo-LTRs showed bursts of formation and recurrent accumulations coinciding with speciation events over past 22 million years. Moreover, our results indicate that the ongoing expansion of ERVK solo-LTRs in these bird species, marked by high transcriptional activity of ERVK retroviral genes in reproductive organs, caused variation of solo-LTRs between individual zebra finches. We experimentally demonstrated that cis-regulatory activity of recently evolved ERVK solo-LTRs may significantly increase the expression level of ITGA2 in the brain of zebra finches compared to chickens. These findings suggest that ERVK solo-LTRs expansion may introduce novel genomic sequences acting as cis-regulatory elements and contribute to adaptive evolution. Overall, our results underscore that the residual sequences of ancient retroviruses could influence the adaptive diversification of species by regulating host gene expression.

Beyond pathogens: the intriguing genetic legacy of endogenous retroviruses in host physiology.

The notion that viruses played a crucial role in the evolution of life is not a new concept. However, more recent insights suggest that this perception might be even more expansive, highlighting the ongoing impact of viruses on host evolution. Endogenous retroviruses (ERVs) are considered genomic remnants of ancient viral infections acquired throughout vertebrate evolution. Their exogenous counterparts once infected the host's germline cells, eventually leading to the permanent endogenization of their respective proviruses. The success of ERV colonization is evident so that it constitutes 8% of the human genome. Emerging genomic studies indicate that endogenous retroviruses are not merely remnants of past infections but rather play a corollary role, despite not fully understood, in host genetic regulation. This review presents some evidence supporting the crucial role of endogenous retroviruses in regulating host genetics. We explore the involvement of human ERVs (HERVs) in key physiological processes, from their precise and orchestrated activities during cellular differentiation and pluripotency to their contributions to aging and cellular senescence. Additionally, we discuss the costs associated with hosting a substantial amount of preserved viral genetic material.

Identification and characterization of endogenous retroviruses upon SARS-CoV-2 infection.

Endogenous retroviruses (ERVs) derived from the long terminal repeat (LTR) family of transposons constitute a significant portion of the mammalian genome, with origins tracing back to ancient viral infections. Despite comprising approximately 8% of the human genome, the specific role of ERVs in the pathogenesis of COVID-19 remains unclear. In this study, we conducted a genome-wide identification of ERVs in human peripheral blood mononuclear cells (hPBMCs) and primary lung epithelial cells from monkeys and mice, both infected and uninfected with SARS-CoV-2. We identified 405, 283, and 206 significantly up-regulated transposable elements (TEs) in hPBMCs, monkeys, and mice, respectively. This included 254, 119, 68, and 28 ERVs found in hPBMCs from severe and mild COVID-19 patients, monkeys, and transgenic mice expressing the human ACE2 receptor (hACE2) and infected with SARS-CoV-2. Furthermore, analysis using the Genomic Regions Enrichment of Annotations Tool (GREAT) revealed certain parental genomic sequences of these up-regulated ERVs in COVID-19 patients may be involved in various biological processes, including histone modification and viral replication. Of particular interest, we identified 210 ERVs specifically up-regulated in the severe COVID-19 group. The genes associated with these differentially expressed ERVs were enriched in processes such as immune response activation and histone modification. HERV1_I-int: ERV1:LTR and LTR7Y: ERV1:LTR were highlighted as potential biomarkers for evaluating the severity of COVID-19. Additionally, validation of our findings using RT-qPCR in Bone Marrow-Derived Macrophages (BMDMs) from mice infected by HSV-1 and VSV provided further support to our results. This study offers insights into the expression patterns and potential roles of ERVs following viral infection, providing a valuable resource for future studies on ERVs and their interaction with SARS-CoV-2.

Roles of Human Endogenous Retrovirus-K-Encoded Np9 in Human Diseases: A Small Protein with Big Functions.

Human Endogenous Retrovirus Sequences (HERVs) constitute up to 8% of the human genome, yet not all HERVs remain silent passengers within our genomes. Some HERVs, especially HERV type K (HERV-K), have been found to be frequently transactivated in a variety of inflammatory diseases and human cancers. Np9, a small protein translated from the HERV-K env reading frame, has been reported as an oncogenic protein and is present in a variety of tumors and transformed cells. The Np9 protein can crosstalk with many cellular factors and is involved in the pathogenicity of various diseases, including some oncogenic virus infections. In the current review, we summarize recent findings about Np9 clinical relevance/implications, its mediated cellular functions/mechanisms, and potential targeted therapies in development.

Human endogenous retrovirus type K encoded Np9 oncoprotein induces DNA damage response.

Human endogenous retrovirus sequences (HERVs) constitute up to 8% of the human genome, yet not all HERVs remain silent passengers within our genomes. Some HERVs, especially the HERV type K (HERV-K), have been found to be frequently transactivated in a variety of inflammatory diseases and human cancers. Np9, a 9-kDa HERV-K encoded protein, has been reported as an oncoprotein and found present in a variety of tumors and transformed cells. In the current study, we for the first time reported that ectopic expression of Np9 protein was able to induce DNA damage response from host cells especially through upregulation of γH2AX. Furthermore, we found that direct knockdown of Np9 by RNAi in Kaposi's Sarcoma-associated herpesvirus (KSHV) infected cells effectively reduced LANA expression, the viral major latent oncoprotein in vitro and in vivo, which may represent a novel strategy against virus-associated malignancies.

Contrasting segregation patterns among endogenous retroviruses across the koala population.

Koalas (Phascolarctos cinereus) have experienced a history of retroviral epidemics leaving their trace as heritable endogenous retroviruses (ERVs) in their genomes. A recently identified ERV lineage, named phaCin-ß, shows a pattern of recent, possibly current, activity with high insertional polymorphism in the population. Here, we investigate geographic patterns of three focal ERV lineages of increasing estimated ages, from the koala retrovirus (KoRV) to phaCin-ß and to phaCin-ß-like, using the whole-genome sequencing of 430 koalas from the Koala Genome Survey. Thousands of ERV loci were found across the population, with contrasting patterns of polymorphism. Northern individuals had thousands of KoRV integrations and hundreds of phaCin-ß ERVs. In contrast, southern individuals had higher phaCin-ß frequencies, possibly reflecting more recent activity and a founder effect. Overall, our findings suggest high ERV burden in koalas, reflecting historic retrovirus-host interactions. Importantly, the ERV catalogue supplies improved markers for conservation genetics in this endangered species.

Activation of endogenous retroviruses and induction of viral mimicry by MEK1/2 inhibition in pancreatic cancer.

While pancreatic ductal adenocarcinomas (PDACs) are addicted to KRAS-activating mutations, inhibitors of downstream KRAS effectors, such as the MEK1/2 kinase inhibitor trametinib, are devoid of therapeutic effects. However, the extensive rewiring of regulatory circuits driven by the attenuation of the KRAS pathway may induce vulnerabilities of therapeutic relevance. An in-depth molecular analysis of the transcriptional and epigenomic alterations occurring in PDAC cells in the initial hours after MEK1/2 inhibition by trametinib unveiled the induction of endogenous retroviruses (ERVs) escaping epigenetic silencing, leading to the production of double-stranded RNAs and the increased expression of interferon (IFN) genes. We tracked ERV activation to the early induction of the transcription factor ELF3, which extensively bound and activated nonsilenced retroelements and synergized with IRF1 (interferon regulatory factor 1) in the activation of IFNs and IFN-stimulated genes. Trametinib-induced viral mimicry in PDAC may be exploited in the rational design of combination therapies in immuno-oncology.

Species-Specific Transcription Factors Associated with Long Terminal Repeat Promoters of Endogenous Retroviruses: A Comprehensive Review.

Endogenous retroviruses (ERVs) became a part of the eukaryotic genome through endogenization millions of years ago. Moreover, they have lost their innate capability of virulence or replication. Nevertheless, in eukaryotic cells, they actively engage in various activities that may be advantageous or disadvantageous to the cells. The mechanisms by which transcription is triggered and implicated in cellular processes are complex. Owing to the diversity in the expression of transcription factors (TFs) in cells and the TF-binding motifs of viruses, the comprehensibility of ERV initiation and its impact on cellular functions are unclear. Currently, several factors are known to be related to their initiation. TFs that bind to the viral long-terminal repeat (LTR) are critical initiators. This review discusses the TFs shown to actively associate with ERV stimulation across species such as humans, mice, pigs, monkeys, zebrafish, Drosophila, and yeast. A comprehensive summary of the expression of previously reported TFs may aid in identifying similarities between animal species and endogenous viruses. Moreover, an in-depth understanding of ERV expression will assist in elucidating their physiological roles in eukaryotic cell development and in clarifying their relationship with endogenous retrovirus-associated diseases.

Correlation analysis of cancer stem cell marker CD133 and human endogenous retrovirus (HERV)-K env in SKOV3 ovarian cancer cells.

BACKGROUND: Human endogenous retrovirus (HERV)-K is a type of retrovirus that is present in the human genome, and its expression is usually silenced in healthy tissues. The precise mechanism by which HERV-K env influences cancer stemness is not fully understood, but it has been suggested that HERV-K env may activate various signaling pathways that promote stemness traits in cancer cells. OBJECTIVE: To establish the connection between HERV-K env expression and cancer stemness in ovarian cancer cells, we carried out correlation analyses between HERV-K env and the cancer stem cell (CSC) marker known as the cluster of differentiation 133 (CD133) gene in SKOV3 ovarian cancer cells. METHOD: To perform correlation analysis between HERV-K env and CSCs, ovarian cancer cells were cultured in a medium designed for cancer stem cell induction. The expression of HERV-K env and CD133 genes was verified using quantitative real-time polymerase chain reaction (RT-qPCR) and Western blot analyses. Additionally, the expression of stemness-related markers, such as OCT-4 and Nanog, was also confirmed using RT-qPCR. RESULTS: In the stem cell induction medium, the number of tumorsphere-type SKOV3 cells increased, and the expression of CD133 and HERV-K env genes was up-regulated. Additionally, other stemness-related markers like OCT-4 and Nanog also exhibited increased expression when cultured in the cancer stem cell induction medium. However, when HERV-K env knockout (KO) SKOV3 cells were cultured in the same cancer stem cell induction medium, there was a significant decrease in the number of tumorsphere-type cells compared to mock SKOV3 cells subjected to the same conditions. Furthermore, the expression of CD133, Nanog, and OCT-4 did not show a significant increase in HERV-K env KO SKOV3 cells compared to mock SKOV3 cells cultured in the same cancer stem cell induction medium. CONCLUSION: These findings indicate that the expression of HERV-K env increased in SKOV3 cells when cultured in cancer stem cell induction media, and cancer stem cell induction was inhibited by KO of HERV-K env in SKOV3 cells. These results suggest a strong association between HERV-K env and stemness in SKOV3 ovarian cancer cells.

Different Host-Endogenous Retrovirus Relationships between Mammals and Birds Reflected in Genome-Wide Evolutionary Interaction Patterns.

Mammals and birds differ largely in their average endogenous retrovirus loads, namely the proportion of endogenous retrovirus in the genome. The host-endogenous retrovirus relationships, including conflict and co-option, have been hypothesized among the causes of this difference. However, there has not been studies about the genomic evolutionary signal of constant host-endogenous retrovirus interactions in a long-term scale and how such interactions could lead to the endogenous retrovirus load difference. Through a phylogeny-controlled correlation analysis on ∼5,000 genes between the dN/dS ratio of each gene and the load of endogenous retrovirus in 12 mammals and 21 birds, separately, we detected genes that may have evolved in association with endogenous retrovirus loads. Birds have a higher proportion of genes with strong correlation between dN/dS and the endogenous retrovirus load than mammals. Strong evidence of association is found between the dN/dS of the coding gene for leucine-rich repeat-containing protein 23 and endogenous retrovirus load in birds. Gene set enrichment analysis shows that gene silencing rather than immunity and DNA recombination may have a larger contribution to the association between dN/dS and the endogenous retrovirus load for both mammals and birds. The above results together showing different evolutionary patterns between bird and mammal genes can partially explain the apparently lower endogenous retrovirus loads of birds, while gene silencing may be a universal mechanism that plays a remarkable role in the evolutionary interaction between the host and endogenous retrovirus. In summary, our study presents signals that the host genes might have driven or responded to endogenous retrovirus load changes in long-term evolution.

Development of human endogenous retrovirus type K- related treatments for human diseases.

Human endogenous retroviruses (HERVs) constitute approximately 8% of the human genome and have long been regarded as silent passengers within our genomes. However, the reactivation of HERVs has been increasingly linked to a range of human diseases, particularly the HERV-K (HML-2) family. Many studies are dedicated to elucidating the potential role of HERV-K in pathogenicity. While the underlying mechanisms require further investigation, targeting HERV-K transactivation emerges as a promising avenue for treating human diseases, including cancer, autoimmune disorders, neurodegenerative conditions, and infectious diseases. In this review, we summarize recent advancements in the development of HERV-K-targeted therapeutic strategies against various human diseases, including antiretroviral drugs, immunotherapy, and vaccines.

Evolution of Endogenous Retroviruses in the Subfamily of Caprinae.

The interest in endogenous retroviruses (ERVs) has been fueled by their impact on the evolution of the host genome. In this study, we used multiple pipelines to conduct a de novo exploration and annotation of ERVs in 13 species of the Caprinae subfamily. Through analyses of sequence identity, structural organization, and phylogeny, we defined 28 ERV groups within Caprinae, including 19 gamma retrovirus groups and 9 beta retrovirus groups. Notably, we identified four recent and potentially active groups prevalent in the Caprinae genomes. Additionally, our investigation revealed that most long noncoding genes (lncRNA) and protein-coding genes (PC) contain ERV-derived sequences. Specifically, we observed that ERV-derived sequences were present in approximately 75% of protein-coding genes and 81% of lncRNA genes in sheep. Similarly, in goats, ERV-derived sequences were found in approximately 74% of protein-coding genes and 75% of lncRNA genes. Our findings lead to the conclusion that the majority of ERVs in the Caprinae genomes can be categorized as fossils, representing remnants of past retroviral infections that have become permanently integrated into the genomes. Nevertheless, the identification of the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 groups indicates the presence of relatively recent and potentially active ERVs in these genomes. These particular groups may contribute to the ongoing evolution of the Caprinae genome. The identification of putatively active ERVs in the Caprinae genomes raises the possibility of harnessing them for future genetic marker development.

Endogenous retrovirus HERVH-derived lncRNA UCA1 controls human trophoblast development.

Endogenous retroviruses (ERVs) are frequently reactivated in mammalian placenta. It has been proposed that ERVs contribute to shaping the gene regulatory network of mammalian trophoblasts, dominantly acting as species- and placental-specific enhancers. However, whether and how ERVs control human trophoblast development through alternative pathways remains poorly understood. Besides the well-recognized function of human endogenous retrovirus-H (HERVH) in maintaining pluripotency of early human epiblast, here we present a unique role of HERVH on trophoblast lineage development. We found that the LTR7C/HERVH subfamily exhibits an accessible chromatin state in the human trophoblast lineage. Particularly, the LTR7C/HERVH-derived Urothelial Cancer Associated 1 (UCA1), a primate-specific long non-coding RNA (lncRNA), is transcribed in human trophoblasts and promotes the proliferation of human trophoblast stem cells (hTSCs), whereas its ectopic expression compromises human trophoblast syncytialization coinciding with increased interferon signaling pathway. Importantly, UCA1 upregulation is detectable in placental samples from early-onset preeclampsia (EO-PE) patients and the transcriptome of EO-PE placenta exhibits considerable similarities to that of the syncytiotrophoblasts differentiated from UCA1-overexpressing hTSCs, supporting up-regulated UCA1 as a potential biomarker of this disease. Altogether, our data shed light on the versatile regulatory role of HERVH in early human development and provide a unique mechanism whereby ERVs exert a function in human placentation and placental syndromes.

The expression analysis of human endogenous retrovirus-K Env, Np9, and Rec transcripts in cervical cancer.

While infection with high-risk human papillomavirus (HPV) types is necessary for cervical cancer (CC) development, it is not enough, and other risk factors are required. Several studies have reported the activation of HERV-K in different cancers; however, the investigation of HERV-K expression levels in CC is scarce. In this study, it was hypothesized that activation of HERV-K could play an essential role in CC development. In this order, the expression levels of HERV-K Env, Np9, and Rec transcripts were investigated on 147 normal to CC uterine cervical tissues using quantitative real-time PCR. The significantly higher levels of HERV-K Env and Np9 transcripts were found in patients with cervical intraepithelial neoplasia (CIN) II-III and CC groups compared to those in the normal/CIN I group. Expression of Rec transcript was also higher only in the CC group than normal/CIN I group. Among CC patients, meaningfully higher levels of HERV-K Env and Np9 transcripts were found in patients with squamous cell carcinoma rather than in adenocarcinoma. When only the HPV 16 positive samples were investigated, it was found that the mean difference in Env and Np9 mRNA levels was meaningfully higher among precancer lesions and the cancer group in comparison with the normal group. However, the Rec mRNA level showed no significant differences. The association between the expression of HERV-K genes was investigated, and a significant positive correlation of Env expression with Np9 transcript was found only in the group with precancer lesions (R = 0.6, p = 0.0037). Moreover, a significant positive correlation was found between Rec and Np9 transcripts in patients with normal cervix tissues (R = 0.26, p = 0.033). However, no correlations were observed between the expression of Env and Rec in the three groups. In conclusion, our results showed that HERV-K transcripts, especially Env and Np9, upregulated during cervical lesion progression. These findings highlight the potential use of HERV-K Env and Np9 as biomarkers for CC diagnosis and prognosis. Further investigation is needed to determine the clinical utility of these markers and whether targeting HERV-K oncogenes could be a viable therapeutic strategy for CC.

GWAS reveals determinants of mobilization rate and dynamics of an active endogenous retrovirus of cattle.

Five to ten percent of mammalian genomes is occupied by multiple clades of endogenous retroviruses (ERVs), that may count thousands of members. New ERV clades arise by retroviral infection of the germline followed by expansion by reinfection and/or retrotransposition. ERV mobilization is a source of deleterious variation, driving the emergence of ERV silencing mechanisms, leaving "DNA fossils". Here we show that the ERVK[2-1-LTR] clade is still active in the bovine and a source of disease-causing alleles. We develop a method to measure the rate of ERVK[2-1-LTR] mobilization, finding an average of 1 per ~150 sperm cells, with >10-fold difference between animals. We perform a genome-wide association study and identify eight loci affecting ERVK[2-1-LTR] mobilization. We provide evidence that polymorphic ERVK[2-1-LTR] elements in four of these loci cause the association. We generate a catalogue of full length ERVK[2-1-LTR] elements, and show that it comprises 15% of C-type autonomous elements, and 85% of D-type non-autonomous elements lacking functional genes. We show that >25% of the variance of mobilization rate is determined by the number of C-type elements, yet that de novo insertions are dominated by D-type elements. We propose that D-type elements act as parasite-of-parasite gene drives that may contribute to the observed demise of ERV elements.

Comprehensive characterization of ERV-K (HML-8) in the chimpanzee genome revealed less genomic activity than humans.

Endogenous retroviruses (ERVs) originate from ancestral germline infections caused by exogenous retroviruses. Throughout evolution, they have become fixed within the genome of the animals into which they were integrated. As ERV elements coevolve with the host, they are normally epigenetically silenced and can become upregulated in a series of physiological and pathological processes. Generally, a detailed ERV profile in the host genome is critical for understanding the evolutionary history and functional performance of the host genome. We previously characterized and cataloged all the ERV-K subtype HML-8 loci in the human genome; however, this has not been done for the chimpanzee, the nearest living relative of humans. In this study, we aimed to catalog and characterize the integration of HML-8 in the chimpanzee genome and compare it with the integration of HML-8 in the human genome. We analyzed the integration of HML-8 and found that HML-8 pervasively invaded the chimpanzee genome. A total of 76 proviral elements were characterized on 23/24 chromosomes, including detailed elements distribution, structure, phylogeny, integration time, and their potential to regulate adjacent genes. The incomplete structure of HML-8 proviral LTRs will undoubtedly affect their activity. Moreover, the results indicated that HML-8 integration occurred before the divergence between humans and chimpanzees. Furthermore, chimpanzees include more HML-8 proviral elements (76 vs. 40) and fewer solo long terminal repeats (LTR) (0 vs. 5) than humans. These results suggested that chimpanzee genome activity is less than the human genome and that humans may have a better ability to shape and screen integrated proviral elements. Our work is informative in both an evolutionary and a functional context for ERVs.

Genome-wide characterization and comparison of endogenous retroviruses among 3 duck reference genomes.

Endogenous retroviruses (ERV) are viral genomes integrated into the host genome and can be stably inherited. Although ERV sequences have been reported in some avian species' genome, the duck endogenous retroviruses (DERV) genome has yet to be quantified. This study aimed to identify ERV sequences and characterize genes near ERVs in the duck genome by utilizing LTRhavest and LTRdigest tools to forecast the duck genome and analyze the distribution of ERV copies. The results revealed 1,607, 2,031, and 1,908 full-length ERV copies in the Pekin duck (ZJU1.0), Mallard (CAU_wild_1.0), and Shaoxing duck (CAU_laying_1.0) genomes, respectively, with average lengths of 7,046, 7,027, and 6,945 bp. ERVs are mainly distributed on the 1, 2, and sex chromosomes. Phylogenetic analysis demonstrated the presence of Betaretrovirus in 3 duck genomes, whereas Alpharetrovirus was exclusively identified in the Shaoxing duck genome. Through screening, 596, 315, and 343 genes adjacent to ERV were identified in 3 duck genomes, respectively, and their functions of ERV neighboring genes were predicted. Functional enrichment analysis of ERV-adjacent genes revealed enrichment for Focal adhesion, Calcium signaling pathway, and Adherens junction in 3 duck genomes. The overlapped genes were highly expressed in 8 tissues (brain, fat, heart, kidney, liver, lung, skin, and spleen) of 8-wk-old Mallard, revealing their important expression in different tissues. Our study provides a new perspective for understanding the quantity and function of DERVs, and may also provide important clues for regulating nearby genes and affecting the traits of organisms.